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Introduction  and  Background 

Many  turbine  engine  components  are  subjected  to  complex  loadings  in  service  that  vary  in 
direction  and  intensity.  Accurate  and  reliable  life  predictions  for  such  components  require 
methodologies  that  can  assess  the  damage  caused  by  individual  load  events  and  provide  a  means 
for  estimating  the  accumulation  of  damage  under  realistic  loading  scenarios.  The  consideration 
of  important  contributing  factors,  including  multiaxiality  of  the  stress  field,  interactions  between 
LCF  and  HCF  loads,  and  the  rate  of  damage  accumulation  under  variable  loading  conditions, 
presents  challenges  in  the  development  of  a  robust,  yet  useable  design  methodology. 
Consequently,  most  component  and  system  lifmg  procedures  currently  in  use  have  adopted 
simplifying  assumptions,  coupled  with  experimental  verification  through  elaborate  component 
and  engine  testing  regiments,  to  establish  durability  limits  for  aircraft  engines. 

Despite  a  substantial  body  of  research  showing  that  the  accumulation  of  fatigue  damage  in 
metals  is  dependent  on  several  factors,  including  load  sequence,  load  path,  and  load  interactions, 
the  Palmgren-Miner  linear  damage  rule  (Miner’s  rule)  is  still  used  almost  exclusively  in 
commercial  fatigue  life  prediction  algorithms.  Although  it  is  generally  recognized  that  the  linear 
damage  assumption  may  lead  to  inaccuracies,  it  is  utilized  for  its  computational  simplicity  and 
due  to  the  widespread  recognition  that  there  are  no  other  proven  methodologies  capable  of 
providing  consistently  better  results.  Numerous  nonlinear  cumulative  damage  models  have  been 
proposed  in  the  technical  literature.  However,  these  models  have  generally  been  derived  in 
response  to  specific  phenomena  observed  under  limited  experimental  test  programs,  most 
notably  under  uniaxial  loading  conditions.  Such  models,  while  often  showing  good  correlation 
to  a  specific  data  set,  have  not  demonstrated  any  widespread  ability  to  model  fatigue  damage 
accumulation  under  a  variety  of  loading  conditions.  These  models  are,  for  the  most  part, 
empirical  approaches  that  have  relied  little  on  analytical  considerations  of  the  causes  of  nonlinear 
damage  accumulation. 

Previous  studies  performed  by  the  PI  [1  -  3],  demonstrated  that  fatigue  damage  accumulation 
in  aircraft  engine  materials  is  highly  sensitive  to  load  path  (multiaxiality),  load-level  (LCF/HCF) 
and  load-path  interactions  among  cycles  in  the  history,  and  the  relative  amounts  of  LCF  damage 
and  HCF  damage  in  the  load  history.  Thus,  a  robust  lifing  methodology  for  turbine  engine 
components  must  incorporate  the  following: 
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•  A  multiaxial  fatigue  damage  parameter  capable  of  accurately  predicting  fatigue  lives 
under  a  variety  of  loading  conditions,  including  non-proportional  loadings  and  multiaxial 
mean  stresses,  over  a  range  of  levels  spanning  LCF  to  HCF  loadings.  This  parameter 
must  be  developed  based  on  a  physical  interpretation  of  the  damage  mechanisms  so  that  it 
possesses  the  capability  to  distinguish  between  the  different  damage  mechanisms 
responsible  for  the  nucleation  and  growth  of  small  cracks  in  the  LCF  and  HCF  regimes. 

•  A  means  for  defining  and  identifying  individual  load  cycles  within  a  multiaxial,  variable 
load  sequence.  While  accepted  cycle-counting  procedures  have  long  been  recognized  for 
uniaxial  load  sequences,  the  extension  to  multiaxial  loadings  brings  forth  additional 
challenges,  particularly  when  stress  components  cycle  out  of  phase  of  one  another. 

•  A  cumulative  damage  rule  capable  of  accounting  for  load-sequence  effects,  load 
interactions  among  LCF  and  HCF  cycles,  and  nonlinear  damage  accumulation  rates  that 
are  dependent  on  the  degree  of  LCF  and  HCF  damage  in  the  service  history. 

•  Finally,  the  methodology  must  be  easily  implemented  and  adaptable  to  current  lifing 
algorithms  that  are  utilized  by  design  engineers. 

This  report  summarizes  the  results  of  a  study  focusing  on  the  development  of  a  robust 
methodology  for  cumulative  fatigue  damage  analysis  of  aircraft  engine  alloys.  This  effort 
consisted  of  an  experimental  component  and  an  analytical  component.  The  experimental  test 
program  was  designed  to  provide  critical  information  related  to  damage  accumulation,  including 
effects  such  as  multiaxiality  of  the  stress  field  and  LCF/HCF  interactions.  A  new  critical-plane 
based  multiaxial  fatigue  damage  parameter  was  proposed  that  provided  good  correlation  with 
experimental  test  data.  Finally,  a  nonlinear  cumulative  damage  model  was  developed,  based  on 
the  damage  curve  concept,  that  accounts  for  load-sequence  and  load-interaction  effects,  as  well 
as  damage  caused  by  sub-threshold  HCF  cycles.  The  details  of  the  experimental  and  analytical 
results  are  provided  in  the  following  sections. 

Experimental  Results:  Multiaxial  Fatigue  Tests 

A  total  of  62  multiaxial  fatigue  specimens  were  machined  from  forged  plates  of  Ti-6A1-4V. 
This  material  was  chosen  for  analysis  because  it  is  commonly  used  in  aircraft  turbine  engines,  is 
readily  available,  and  has  been  well  characterized.  Both  tubular  (hollow)  and  solid  specimens 
were  utilized  in  this  study.  Thin-walled  tubes  are  preferred  for  tension/torsion  tests  due  to  the 
simplified  analysis  that  results  from  the  assumption  of  constant  shear  stress  through  the  tube 
walls.  Provided  the  applied  loads  are  not  too  large,  as  is  typical  of  HCF  loading  conditions,  thin- 
walled  tubes  can  be  safely  tested  in  tension/torsion.  However,  when  subjected  to  higher  loads 
representative  of  LCF  conditions,  there  is  a  concern  of  thin-walled  buckling  in  torsion. 
Consequently,  several  biaxial  LCF  tests  were  conducted  using  solid  specimens. 

The  test  section  for  the  tubular  specimens  had  a  nominal  OD  of  0.375  in  and  ID  of  0.25  in, 
and  a  length  of  1.5  in.  The  solid  specimens  had  a  nominal  test  section  diameter  of  0.5  in  and 
length  of  1.2  in.  All  specimens  were  subjected  to  a  low  stress  grind  to  16  RMS,  and  then 
longitudinally  polished  to  8  RMS.  A  total  of  52  tubular  specimens  and  12  solid  specimens  were 
machined. 

All  multiaxial  fatigue  tests  were  performed  at  the  Structural  Integrity  Laboratory  at  the 
University  of  Utah,  under  the  direction  of  Dr.  Dan  Adams  and  Jeff  Kessler,  and  the  Advanced 
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Materials  Testing  Laboratory  at  the  University  of  Illinois,  under  the  direction  of  Dr.  Peter 
Kurath.  The  LCF  tests  were  conducted  in  strain  control  due  to  the  high  level  of  yielding  along 
the  outer  surface  of  the  specimens.  The  HCF  tests,  which  experienced  very  low  levels  of  plastic 
behavior,  were  conducted  under  load-control  conditions.  An  elastic-plastic  finite  element 
analysis  was  performed  to  determine  surface  stresses  in  specimens  that  incurred  plastic  strains  at 
the  peak  load  level. 

Several  constant-amplitude  biaxial  fatigue  tests  were  first  performed  to  assist  in  the 
development  of  the  multiaxial  fatigue  damage  parameter,  as  well  as  to  calibrate  the  final  model. 
These  tests  consisted  of  uniaxial,  torsion,  and  proportional  tension/torsion  tests  at  stress  ratios  of 
R  =  -1  and  R  =  0.  Load  levels  were  varied  to  produce  fatigue  lives  covering  the  spectrum  from 
LCF  to  HCF  conditions.  The  results  of  the  constant-amplitude  tests  are  listed  in  Appendix  1. 
The  elastic-plastic  stresses  listed  in  this  table  represent  the  surface  stresses  at  maximum  and 
minimum  loads  as  determined  by  an  elastic-plastic  finite  element  analysis  using  ANSYS. 

Next,  a  series  of  simulated  “mission”  tests  were  conducted,  consisting  of  a  mix  of  LCF  and 
HCF  cycles  with  varying  load  paths.  These  tests  were  conducted  to  analyze  the  effects  of  load- 
path  and  load-level  interactions.  The  load  paths  used  in  this  phase  of  the  test  program  were  the 
same  as  those  used  in  the  constant-amplitude  tests.  However,  the  types  of  load  paths  were  varied 
between  LCF  and  HCF  cycles  in  order  to  create  differing  relationships  in  damage  mechanisms 
among  the  cycles  in  these  tests.  The  intent  was  to  better  understand  how  variations  in  load  path 
between  cycles  affect  the  rate  of  damage  accumulation  under  multiaxial  loading  conditions. 

All  multiaxial  mission  tests  were  comprised  of  a  block  of  cycles  consisting  of  1  LCF  cycle 
coupled  with  50  HCF  cycles,  as  illustrated  schematically  for  uniaxial  loading  conditions  in  Fig. 
1.  In  all  cases,  LCF  lives  were  on  the  order  of  10?  cycles,  while  HCF  levels  were  at  or  below 
threshold  levels  (infinite  life  conditions).  Thus,  on  the  basis  of  a  conventional  linear  cumulative 
damage  analysis,  the  HCF  cycles  would  contribute  no  damage,  meaning  the  mission  lives  would 
be  expected  to  be  similar  to  LCF  lives.  The  load  paths  used  in  these  tests  are  illustrated 
schematically  in  Fig.  2.  The  test  conditions  and  results  are  listed  in  Appendix  2. 


Figure  1. 
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Schematic  illustration  of  mission  histories. 
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Figure  2.  Multiaxial  load  paths  used  in  the  mission  histories. 

As  previously  noted,  stress  and  strain  values  on  the  surface  of  the  specimens  experiencing 
yielding  were  determined  from  an  elastic-plastic  finite  element  analysis  using  ANSYS.  This 
analysis  employed  a  multilinear  kinematic  hardening  rule  in  conjunction  with  the  Ramberg- 
Osgood  cyclic  stress-strain  curve,  as  measured  from  strain-controlled  uniaxial  tests.  This  curve 
is  shown  in  Fig.  3,  along  with  the  experimental  data.  The  material  parameters  for  the  Ramberg- 
Osgood  formulation,  shown  in  Eq.  (1),  were  found  to  be  E  =  111.8  GPa,  K’  =  1383  MPa,  and 
n’  =  0.0956. 
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Figure  3.  Cyclic  stress-strain  curve. 

The  results  of  the  mission  tests  are  presented  graphically  in  Fig.  4,  relative  to  the 
corresponding  LCF  tests.  This  figure  plots  the  mission  lives  (square  symbols)  in  relation  to  the 
LCF  lives  (circles)  for  a  variety  of  load  paths  (tensile,  torsion,  and  proportional).  As  is  evident 
from  these  results,  the  presence  of  very  small  HCF  cycles  in  a  load  history  caused  a  substantial 
reduction  in  the  mission  lives,  relative  to  LCF  lives,  in  some  cases.  Notably,  some  mission  lives 
were  reduced  by  a  factor  of  20  or  more  due  to  the  HCF  cycles;  however,  the  life  reductions  were 
heavily  dependent  on  the  load  path  relationships  between  the  LCF  and  HCF  cycles. 
Conventional  damage  accumulation  techniques,  such  as  those  based  on  the  Palmgren-Miner 
linear  damage  rule,  would  thus  drastically  overestimate  the  mission  lives.  These  results  indicate 
a  pressing  need  to  develop  more  accurate  life  prediction  techniques  for  variable  amplitude 
loading  that  take  into  account  nonlinear  damage  accumulation  and  load  path  interactions. 
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Figure  4.  Comparison  of  mission  lives  and  LCF  lives. 


Multiaxial  Fatigue  Modeling:  Critical  Plane  Analysis 

The  approach  taken  in  this  study  for  cumulative  damage  analysis  of  aircraft  engine 
components  focuses  on  the  nucleation  and  growth  of  small  cracks  (microcracks),  rather  than  the 
propagation  of  macrocracks.  Many  engine  components  are  currently  life-limited  on  the  basis  of 
crack  nucleation,  as  this  typically  consumes  the  majority  of  the  total  life  under  the  loading 
conditions  to  which  these  components  are  subjected.  Since  a  requirement  of  the  damage 
accumulation  model  is  that  it  be  easy  to  use  by  a  design  engineer,  it  is  important  that  this  model 
be  capable  of  being  integrated  within  existing  design  algorithms.  Furthermore,  the  effect  of  load 
multiaxiality  on  crack  growth  appears  to  be  less  of  a  concern  than  on  crack  nucleation.  Many 
researchers  have  reported  that,  under  mixed-mode  loading,  macrocracks  tend  to  grow 
predominantly  under  mode-I  conditions.  In  addition,  conventional  analysis  techniques  for  crack 
growth,  based  on  fracture  mechanics  principles,  inherently  account  for  nonlinearity  in  damage 
growth  and  load-sequence  effects  through  an  incremental  crack  propagation  analysis  that  can 
model  crack  retardation  due  to  overloads  using  crack  closure  and  plastic  zone  concepts. 
Conversely,  for  crack  nucleation,  there  are  no  broadly  applicable  techniques  currently  available 
for  cumulative  damage  analysis  other  than  Miner’s  rule,  despite  overwhelming  experimental 
evidence  of  nonlinear  damage  accumulation.  Thus,  the  focus  of  this  study  was  to  provide  a 
mission-based  lifing  methodology  for  engine  components  based  on  crack  nucleation  concepts 
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that  accounts  for  multiaxial  stress  states  and  nonlinear  damage  accumulation,  and  which  can  be 
implemented  into  existing  design  codes. 

Past  research  by  the  PI  [4]  resulted  in  the  development  of  a  critical-plane  parameter  for  crack 
nucleation  that  demonstrated  an  excellent  ability  to  correlate  uniaxial  and  multiaxial 
(tension/torsion)  fatigue  data  for  aircraft  engine  alloys,  including  Ti-6A1-4V.  This  parameter 
was  expressed  as 


DP  =  t„ 
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In  Eq.  (2),  xmax  and  Tmin  are  the  shear  reversal  points  on  the  critical  plane  (defined  here  as  the 
plane  of  maximum  alternating  shear  stress),  and  a-max  and  Gxmin  are  the  corresponding  normal 
stresses  on  this  plane  at  the  instant  of  shear  reversal.  These  stresses  are  divided  by  the  yield 
strength,  Gy,  to  maintain  consistency  of  units.  In  the  second  term,  Gmax  and  Gmin  are  the 
maximum  and  minimum  normal  stresses  on  the  critical  plane  over  the  entire  cycle.  The  c 
notation  indicates  that  only  positive  (tensile)  values  are  used  in  the  equation.  If  Gmin  or  Gmax  is 
negative,  that  value  is  set  to  zero  in  Eq.  (2). 

The  derivation  of  this  parameter  was  based  on  a  physical  interpretation  of  the  mechanisms 
governing  the  nucleation  and  growth  of  small  fatigue  cracks.  The  first  term  in  Eq.  (2)  represents 
the  primary  mechanism  for  crack  nucleation  in  ductile  metals;  i.e.,  cyclic  shear  stresses.  The 
form  of  the  expression  allows  for  the  possibility  of  a  mean  shear  stress  effect,  which  has  been 
observed  in  several  aircraft  engine  alloys.  The  multiplier  term  (containing  Gxmax  and  Gxmin)  acts 
to  modify  the  shear  component  to  account  for  crack  face  interaction  due  to  the  normal  stresses  at 
the  shear  reversals.  A  tensile  normal  stress  will  hold  the  crack  open,  reducing  crack  face 
interaction  and  friction,  thereby  accelerating  crack  advance  along  the  shear  plane.  Conversely,  a 
compressive  normal  stress  will  increase  crack  face  interaction  and  friction,  inhibiting  the 
advance  of  the  crack.  The  definition  of  the  ki  term  allows  for  a  different  influence  of  tensile  and 
compressive  stresses.  The  last  term  in  Eq.  (2)  accounts  for  the  additional  damage  caused  by 
cyclic  tensile  stresses  on  the  critical  plane.  If  multiple  normal  stress  cycles  exist  within  the  shear 
cycle,  their  contribution  is  included  through  the  summation  in  this  term.  The  values  of  kj,  Wi,  k2 
and  W2  are  fit  by  collapsing  uniaxial  and  multiaxial  fatigue  data. 

A  key  aspect  inherent  in  the  damage  parameter  presented  in  Eq.  (2)  is  the  potential  to 
distinguish  between  different  damage  mechanisms  (e.g.,  shear  vs.  tensile).  Prior  work  by  the  PI 
has  shown  a  change  in  damage  mechanisms  between  LCF  and  HCF  cycles  under  multiaxial  load 
paths,  which  may  contribute  to  the  highly  nonlinear  rates  of  damage  accumulation  observed  in 
this  study.  In  order  to  accurately  model  such  effects,  it  is  necessary  to  interpret  the  damage 
mechanisms  between  cycles  and  the  potential  for  load-cycle  interactions.  The  parameter  in  Eq. 
(2)  separates  fatigue  damage  into  shear  and  tensile  modes;  however,  there  are  drawbacks 
associated  with  this  parameter  that  render  it  unsuitable  for  more  complex  variable  load  histories. 
For  example,  if  multiple  normal  stress  cycles  exist  within  a  single  shear  cycle,  the  summation  on 
the  last  term  may  overestimate  the  damage  associated  with  the  major  cycle.  Hence,  it  is 
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necessary  to  break  this  damage  term  out  independently  and  account  for  it  within  the  cumulative 
damage  model.  In  this  manner,  the  model  will  not  only  allow  for  nonlinear  damage 
accumulation  from  cycle  to  cycle,  but  also  within  an  individual  cycle  due  to  different  damage 
mechanisms.  Furthermore,  the  separation  of  damage  mechanisms  on  the  critical  plane  allows  the 
ability  to  weight  the  importance  of  those  mechanisms  differently  in  accordance  with  observed 
phenomena;  e.g.,  shear  mechanisms  dominate  within  LCF  cycles  while  tensile  mechanisms 
dominate  within  HCF  cycles. 

In  consideration  of  the  issues  presented  above,  along  with  a  critical  examination  of  constant 
amplitude  multiaxial  fatigue  data  generated  in  this  study  and  previous  studies  by  the  PI,  a  new 
critical  plane  parameter  was  developed  for  LCF  loadings.  This  parameter  is  shown  in  Eq.  (3). 

DP  =  (3) 

In  this  equation,  Ay  is  the  shear  strain  range  on  the  critical  plane  (defined  as  the  plane  of 
maximum  cyclic  shear  stress/strain),  Tmax  is  the  maximum  absolute  value  of  shear  stress  on  this 
plane  (which  allows  the  model  to  account  for  a  mean  shear  stress  effect),  (xa)max  is  the  maximum 
product  of  shear  and  normal  stress  on  the  critical  plane,  Go  is  an  arbitrary  reference  stress  value 
(to  maintain  unit  consistency),  G  is  the  shear  modulus,  and  w  and  k  are  material  parameters 
determined  from  fitting  experimental  data.  The  last  term  in  this  equation  is  designed  to  account 
for  the  influence  of  the  normal  stress  on  the  critical  plane.  A  tensile  normal  stress  will  open  the 
crack,  thereby  reducing  crack  face  interaction  and  allowing  the  crack  to  nucleate  more  easily  by 
the  cyclic  shear  stress.  However,  this  process  is  also  dependent  on  the  phase  of  the  stresses;  i.e., 
when  a  tensile  normal  stress  occurs  in-phase  with  a  large  shear  stress,  it  will  have  a  greater 
influence  than  if  it  occurs  out-of-phase  with  the  shear  stress. 

The  new  parameter  shown  in  Eq.  (3)  has  several  advantages  over  the  previous  form  shown  in 
Eq.  (2).  The  inclusion  of  the  cyclic  shear  strain,  Ay,  rather  than  just  stress  terms,  allows  the 
parameter  to  better  model  cyclic  hardening  or  softening  that  often  occurs  under  LCF  loading.  In 
addition,  the  absence  of  the  summation  of  small  cycle  damage  within  the  dominant  shear  cycle 
negates  the  overestimation  of  this  type  of  damage  that  plagued  the  previous  parameter.  Instead, 
the  contribution  of  small  cycle  (HCF)  damage  can  be  assessed  independently  in  the  cumulative 
damage  model.  Finally,  the  reduction  of  material  parameters  (from  four  to  two)  that  must  be  fit 
from  experimental  data  results  in  a  much  simplified  model. 

The  new  critical-plane  parameter  shown  in  Eq.  (3)  was  fit  to  the  constant-amplitude 
multiaxial  fatigue  data  for  Ti-6A1-4V  generated  in  this  study,  presented  in  Fig.  5.  As  shown, 
there  is  a  large  variety  of  test  data,  including  uniaxial,  torsion,  and  proportional  load  paths  tested 
under  load-control  and  strain-control  conditions.  From  this  analysis,  the  following  material 
parameters  were  obtained:  Go  =  689  MPa,  w  =  0.46,  and  k  =  0.8.  Although  there  is  a  significant 
amount  of  scatter  in  the  experimental  fatigue  data,  it  is  evident  that  this  parameter  collapses  the 
multiaxial  data  well,  resulting  in  a  simple  and  effective  critical -plane  model  for  the  purposes  of 
LCF  life  prediction  of  high-strength  alloys  subjected  to  multiaxial  loading  conditions. 
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Cycles  to  Failure 


Figure  5.  New  damage  parameter  applied  to  multiaxial  fatigue  data  for  Ti-6A1-4V. 

It  is  theorized  here  that  the  observed  reduction  in  mission  life  due  to  the  HCF  cycles  is 
caused  by  the  propagation  of  microcracks  initiated  by  the  LCF  cycles.  As  it  is  well  known  that 
crack  propagation  is  driven  primarily  by  cyclic  tensile  stresses,  a  tensile-stress  based  damage 
parameter  will  be  required  to  model  the  additional  damage  caused  by  the  HCF  stresses.  In  this 
study,  a  Walker-type  parameter  was  used  for  this  purpose.  The  Walker  parameter  assumes  that 
cyclic  tensile  stresses  cause  fatigue  damage,  and  accounts  for  a  mean  stress  effect  through  a 
modification  incorporating  the  stress  ratio  and  an  adjustable  material  parameter,  w.  The  form  of 
the  Walker  parameter  utilized  in  this  study  is  shown  in  Eq.  (4). 

(_  \W 

1  _  Shsl)  (4) 

In  order  to  develop  a  functional  relationship  for  modeling  the  damage  caused  by  the  HCF 
cycles,  the  Walker  parameter  shown  in  Eq.  (4)  was  fit  to  the  uniaxial  fatigue  data  generated  in 
this  study.  The  resulting  curve  fit  is  shown  in  Fig.  6  with  w  =  0.35.  As  is  evident,  the  Walker 
parameter  provides  a  good  representation  of  the  uniaxial  (tensile)  fatigue  data,  accounting  for  the 
effect  of  mean  stresses.  From  this  figure,  it  can  be  seen  that  values  of  the  Walker  parameter 
below  approximately  700  MPa  (100  ksi)  would  not  be  expected  to  contribute  significantly  to 
fatigue  damage.  Thus,  a  value  of  700  MPa  was  used  in  subsequent  analyses  as  a  “threshold” 
value  for  fatigue  damage  contributions. 
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Figure  6.  Walker  parameter  applied  to  uniaxial  fatigue  data  for  Ti-6A1-4V. 

Nonlinear  Damage  Accumulation  Model 

Numerous  empirical  theories  have  been  proposed  to  model  nonlinear  fatigue  damage 
accumulation  [5],  Among  these  models,  the  Damage  Curve  Approach  (DCA),  developed  by 
Manson,  Halford,  and  coworkers  [6-7],  has  received  greater  attention  due  to  its  simplicity  in 
application  and  its  ability  to  account  for  several  observed  phenomena.  In  this  model,  the  rate  of 
damage  accumulation  is  dependent  on  the  load  level  (or  equivalently,  life  to  failure  at  that  load 
level)  and  the  life  fraction,  n/Nf.  Thus,  different  damage  curves  govern  the  damage  growth  at 
different  life  levels.  The  rate  of  damage  accumulation  is  assumed  to  increase  as  the  number  of 
applied  cycles  increases,  and  the  degree  of  nonlinearity  is  greater  at  higher  life  levels  (lower 
loads).  Using  this  model,  the  instantaneous  state  of  fatigue  damage  is  calculated  as 


where  a  is  a  material  parameter  and  Nref  is  the  life  level  at  which  damage  accumulates  linearly, 
generally  taken  as  unity.  The  model  is  shown  schematically  in  Fig.  7,  illustrating  how  damage 
would  be  computed  for  a  three-level  block  loading  history  in  which  the  load  levels  are 
decreasing. 
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Figure  7.  Illustration  of  the  Damage  Curve  Approach. 

A  desirable  feature  of  this  model  is  that  it  captures  the  load-sequence  effect  commonly 
observed  in  two-level  block  loading  tests.  In  addition,  since  the  damage  is  expressed  solely  in 
terms  of  the  fatigue  life  and  life  ratio,  it  is  independent  of  the  damage  parameter  used  in  the 
calculation  of  fatigue  life  once  a  has  been  obtained.  However,  a  limitation  of  this  model  is  that 
it  neglects  any  damage  caused  by  cycles  below  threshold  or  endurance  values;  i.e.,  if  Nf  =  oo, 
D  =  0.  Thus,  for  the  tests  conducted  in  this  study  (in  which  the  HCF  cycles  were  below 
threshold  values),  the  model  would  predict  no  reduction  in  mission  life  due  to  the  HCF  cycles. 

Due  to  the  limitations  of  the  DCA  model,  a  new  nonlinear  damage  accumulation  model  was 
developed  through  a  modification  to  the  original  form.  The  primary  objective  was  to  allow  the 
model  to  capture  the  damage  contributions  from  very  small  HCF  cycles,  even  those  below 
threshold  stress  (endurance)  limits.  In  this  analysis,  the  sub-threshold  HCF  cycles  are  assumed 
to  cause  no  damage  unless  prior  damage  has  occurred,  such  as  through  the  application  of  LCF 
cycles.  Thus,  in  the  absence  of  existing  damage,  the  HCF  cycles  will  cause  no  fatigue  damage; 
i.e.,  these  cycles  are  too  small  to  nucleate  fatigue  cracks  by  themselves.  However,  if  fatigue 
cracks  are  nucleated  through  LCF  cycles,  the  HCF  cycles  will  then  contribute  to  the  further 
development  and  growth  of  the  microcracks. 

In  the  model,  this  is  accomplished  by  defining  a  “threshold”  damage  value,  D*,  for  small 
HCF  cycles.  This  threshold  value  may  vary  with  different  HCF  cycle  levels.  As  long  as  the  total 
accumulated  damage  (D)  from  prior  cycles  remains  below  D*,  a  given  HCF  cycle  will  not 
contribute  additional  damage.  However,  once  the  total  accumulated  damage  from  prior  cycles 
exceeds  D*,  the  associated  HCF  cycles  will  start  to  contribute  damage. 

The  modified  DCA  model  is  illustrated  in  Fig.  8.  In  this  figure,  the  blue  lines  represent 
finite-life  damage  curves;  i.e.,  curves  associated  with  LCF  cycles  corresponding  to  a  specific 
fatigue  life.  The  red  lines  represent  damage  curves  associated  with  sub-threshold  HCF  cycles  of 
varying  magnitudes,  resulting  in  different  D*  values  for  each  load  level.  Whereas  the  LCF 
cycles  may  start  to  accumulate  damage  immediately,  the  HCF  cycles  will  not  contribute  damage 
until  the  corresponding  D[h  value  has  been  surpassed.  Smaller  HCF  cycles  have  higher  Dth 
values  than  larger  HCF  cycles,  but  eventually  all  cycles  will  contribute  damage. 
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Figure  8.  Illustration  of  the  modified  Damage  Curve  Approach. 

Mathematically,  the  damage  function  for  the  new  model  is  expressed  as 

D  =  Dth  +  (1  -  Dch)  X  (-)  (6) 

where  the  threshold  damage  value,  Dth,  is  given  in  terms  of  the  threshold  damage  parameter 
value  (e.g.,  endurance  limit),  DPth, 


In  using  these  equations,  if  the  current  value  of  the  damage  parameter  (DP)  exceeds  DPth,  then 
Dth  -  0  and  D  is  calculated  directly  from  Eq.  (6).  However,  if  DP  <  DPth,  then  D&  is  first 
calculated  from  Eq.  (7),  with  D  subsequently  determined  from  Eq.  (6).  The  constants  a  and  (3 
are  considered  material  parameters  that  must  be  fit  from  experimental  data.  They  govern  the 
degree  of  nonlinearity  inherent  in  the  accumulation  of  fatigue  damage.  A  desirable  feature  of 
this  model  is  that  different  damage  parameters  can  be  used  for  different  types  of  cycles  to  better 
model  different  damage  mechanisms  that  may  occur  between  cycles;  e.g.,  LCF  vs.  HCF  damage 
mechanisms. 
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The  modified  DCA  model  was  applied  to  the  results  of  the  mission  histories  tested  in  this 
study.  In  this  analysis,  the  shear-based  damage  parameter  shown  in  Eq.  (3)  was  used  for  the  LCF 
cycles,  while  the  tensile-based  Walker  parameter  shown  in  Eq.  (4)  was  used  for  the  HCF  cycles. 
Based  on  the  results  shown  in  Fig.  6,  a  value  of  700  MPa  was  used  for  the  HCF  DP*  value  in  Eq. 
(7).  Optimal  values  of  a  and  (3  were  determined  for  each  test  by  matching  experimental  and 
predicted  lives,  then  these  values  were  averaged  to  determine  a  single,  optimal  value  of  each 
constant  for  the  Ti-6A1-4V  material.  From  this  process,  the  following  values  were  determined: 
a  =  0.25,  p  =  0.5. 

Several  damage  curves  for  the  Ti-6A1-4V  material  using  the  aforementioned  material 
parameters  are  shown  in  Fig.  9.  As  before,  the  blue  curves  represent  finite-life  (LCF)  damage 
curves,  for  which  the  damage  parameter  is  calculated  from  Eq.  (3).  The  red  curves  represent 
HCF  damage,  for  which  the  damage  parameter  is  calculated  from  Eq.  (4).  In  these  cases,  the  DP 
values  are  less  than  the  threshold  value  of  700  MPa,  requiring  the  use  of  both  Eqs.  (6)  and  (7)  to 
calculate  the  current  damage  level. 


Figure  9.  Damage  curves  for  Ti-6A1-4V.  DP  values  shown  refer  to  sub-threshold 
values  of  the  Walker  parameter. 

Using  the  new  damage  model,  predicted  lives  were  calculated  for  the  mission  histories  tested 
in  this  study,  and  compared  to  the  experimental  lives.  The  results  are  shown  in  Fig.  10  in  the 
form  of  a  predicted  vs.  experimental  life  plot.  The  corresponding  predictions  using  the  linear 
damage  rule  are  also  included  for  comparison.  In  this  plot,  the  solid  lines  on  either  side  of  the 
main  diagonal  represent  a  scatter  band  of  a  factor  of  two  in  the  life  predictions,  while  the  dashed 
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lines  represent  a  factor  of  ten  scatter  band.  It  is  clear  that  in  nearly  all  cases,  the  linear  rule  was 
highly  non-conservative,  with  predicted  lives  often  exceeding  experimental  lives  by  a  factor  of 
1 0  or  more.  In  only  two  cases  did  the  linear  damage  rule  provide  accurate  life  predictions  within 
a  factor  of  two  of  the  experimental  lives.  In  contrast,  the  proposed  nonlinear  damage  model 
resulted  in  all  life  predictions  falling  within  a  factor  of  ten,  and  half  of  the  predictions  falling 
within  a  factor  of  two  of  experimental  lives.  Thus,  the  new  model  presents  a  clear  improvement 
in  estimating  the  fatigue  damage  caused  by  small  cycles,  as  well  as  accounting  for  the  load-path 
interactions  between  LCF  and  HCF  cycles  through  consideration  of  the  varying  stress  states  on 
the  critical  plane. 


Experimental  Life 


Figure  10.  Predicted  vs.  experimental  lives  for  the  mission  tests  using  the  new  damage 
model.  Life  predictions  are  also  shown  using  the  conventional  linear 
damage  rule. 


Summary 

A  comprehensive  study  was  undertaken  to  investigate  and  model  the  accumulation  of  fatigue 
damage  under  multiaxial  loading,  with  an  emphasis  on  LCF/HCF  interactions,  load-path 
dependence  between  large  and  small  cycles,  and  nonlinear  rates  of  damage  accumulation. 
Experimental  testing  revealed  a  very  strong  effect  on  damage  accumulation  from  HCF  cycles, 
even  those  that  would  typically  be  considered  below  threshold  or  endurance  levels.  However, 
the  influence  of  the  HCF  cycles  is  heavily  dependent  on  the  load  path.  Notably,  load  histories  in 
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which  cyclic  tensile  stresses  from  HCF  cycles  are  maximized  on  or  near  LCF  critical  shear 
planes  appear  to  exhibit  significant  damage  interactions.  Also,  there  is  a  substantial  influence 
from  HCF  cycles  in  torsional  loading.  However,  under  uniaxial  loading  conditions,  the  effect  of 
the  HCF  cycles  was  found  to  be  negligible.  This  is  an  important  point,  as  most  fatigue  damage 
studies  have  utilized  primarily  uniaxial  test  data. 

Based  on  the  experimental  evidence,  a  new  approach  to  modeling  fatigue  damage 
accumulation  was  presented.  The  approach  relies  on  a  multiaxial  critical  plane  analysis  to 
identify  the  plane  of  crack  nucleation  from  the  LCF  cycles,  assumed  here  to  be  shear  dominated. 
A  separate  parameter,  based  on  tensile  stresses,  is  used  to  assess  the  damage  caused  by  HCF 
cycles  on  the  LCF  critical  plane.  The  rate  of  damage  accumulation  is  assumed  to  be  nonlinear, 
and  dependent  on  the  life  (load)  level  of  each  cycle.  A  “threshold  damage”  level  was  introduced 
for  HCF  cycles.  When  the  accumulated  damage  is  below  the  threshold  level  for  any  given  cycle, 
that  cycle  is  assumed  to  cause  no  damage.  However,  once  exceeded  (such  as  from  larger  cycles), 
even  very  small  cycles  can  contribute  to  fatigue  damage  accumulation.  The  new  model  was  used 
to  predict  fatigue  lives  of  simulated  mission  histories  tested  in  this  study,  and  the  results  were 
found  to  be  significantly  improved  over  the  traditional  linear  damage  rule.  Hence,  this  new 
approach  provides  a  more  reliable  way  to  assess  fatigue  damage  accumulation  in  metals  than 
most  current  methodologies,  specifically  accounting  for  multiaxial  loadings  and  LCF/HCF 
interactions. 
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Table  1.  Constant  Amplitude  Fatigue  Data 
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Table  2.  Test  Data  from  Multiaxial  Mission  Histories 


Table  2  (cont).  Test  Data  from  Multiaxial  Mission  Histories 


